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value y s = — 1 for the irrelevant exponent. 
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Finite-size scaling concepts are a powerful tool for the determination of critical properties 
at phase transitions especially when coupled to phenomenological renormalisation [Q] 
and conformal invariance H ideas. Here we present results from numerical transfer- 
matrix calculations of the correlation length, and quantities derived therefrom, for site 
percolation on infinite strips with free boundary conditions (FBC). The use of FBC 
allows one to assess surface critical properties, including so-called special and surface 
transitions || when they occur. 

Our transfer-matrix formulation of the percolation problem relies on the direct 
application of connectivity concepts |J, as opposed to taking the s — > 1 limit of the 
s -state Potts model ||, which corresponds to bond percolation (and should then, by 
universality, give the same exponents as for the site problem). While the latter approach 
benefits from being a systematic expansion in terms of Whitney polynomials, it is devised 
for general, continuous s and thus carries a high degree of inherent complexity. As shown 
in earlier work |^|, [7|, || and below, the geometric picture based on cluster connectivity 
allows for a straightforward algorithm to be built, from which a nicely extrapolating 
sequence of finite-size estimates is extracted. 

We use strips of width L < 10 sites, both for square and triangular lattices. This 
is the same maximum width reached with periodic boundary conditions (PBC) ]/J ||, 
though in the present case lower symmetry implies that the matrices are of considerably 
larger dimension than on a cylindrical geometry. First, standard phenomenological 
renormalisation (PR) [@ is performed on the site occupation probability p between 
strips of consecutive widths, from which estimates of the critical concentration p c , the 
temperature-like exponent y p and the surface decay-of-correlations exponent, r] s are 
obtained. Alternative finite-size sequences for rj s are produced by setting p at the exact 
(or extrapolated) p c . We then search for a special transition, by introducing a distinct 
probability p s for site occupation along the strip boundaries. A two-parameter PR 
analysis is carried out, by comparing correlation lengths on three strips of consecutive 
widths ||, [lO], 11 1 . Only one non-trivial fixed point, shown to correspond to an ordinary 
transition, is found upon numerical examination of the recursion relations. No evidence 
is found for the existence of a multicritical point related to a special, or surface, 
transition. 

The exponent that governs the decay of correlations along the surface of a semi- 
infinite plane at criticality, i] s , is related to the correlation length on a strip with FBC 
by: 

2L 



a result from conformal invariance [|T^, [Td[| . Note that for a triangular lattice with FBC 
the strip width L = where N is the number of sites across the strip. 

Corrections to scaling must be dealt with, since e.g. Equation (|I|) is expected to be 
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valid only asymptotically. Throughout this work, extrapolations toward L — > oo have 
been done using the Bulirsch-Stoer (BST) algorithm |TJ, |15|. As extensively discussed 
elsewhere | I5| , whenever the leading correction-to-scaling exponent uo is not known a 
priori BST extrapolations rely on keeping it as a free parameter within an interval 
guessed to be reasonable. Central estimates and error bars are evaluated self-consistently 
by selecting the range of uj for which overall fluctuations are minimised. In the following 
we have allowed 0.45 < u < 2.4 for all quantities, as the most likely candidates in the 
case are u = 1 and 2 (see below). 

We implement standard, one-parameter, PR in the usual way by looking for the 
fixed point p* of the implicit recursion relation: 

L ~ L-l ' {) 
where = _ l/hiA^p) is given in terms of the largest eigenvalue A L (p) of the 

column-to-column transfer matrix ||; p* is thus a finite-size estimate of p c . At the 
fixed point, the temperature-like exponent y p = \jv is evaluated by taking suitable 
derivatives IJ. For consistency, rj s is obtained from Equation (|l|) with £ calculated at 
p*. 



Table 1. Results from one-parameter PR. Uncertainties in last quoted digits are shown 
in parentheses. Extrapolations obtained by BST algorithm with correction-to-scaling 
exponent u> in ranges shown. Expected values are exact, unless otherwise noted. 



Square Triangular 



L 


P* 


y P 


Vs 


P* 


Vv 


Vs 


3 


0.671130 


0.662822 


0.290469 


0.573092 


0.665596 


0.337283 


4 


0.644177 


0.676427 


0.351082 


0.547377 


0.681787 


0.393936 


5 


0.629524 


0.686306 


0.394701 


0.533471 


0.692689 


0.434325 


G 


0.620566 


0.693798 


0.427621 


0.525062 


0.700557 


0.464391 


7 


0.614644 


0.699685 


0.453367 


0.519561 


0.706500 


0.487600 


8 


0.610502 


0.704439 


0.474064 


0.515749 


0.711147 


0.506047 


9 


0.607478 


0.708362 


0.491070 


0.512990 


0.714869 


0.521047 


10 


0.605197 


0.711655 


0.505297 


0.510928 


0.717880 


0.533458 


Expected 


0.592745(2)" 


3/4 


2/3 h 


1/2 


3/4 


2/3 fc 


Extrapolated 


0.5925(5) 


0.750(2) 


0.666(3) 


0.5005(2) 


0.750(2) 


0.676(3) 




1.50(50) 


1.10(10) 


1.00(5) 


2.00(10) 


1.10(10) 


1.05(5) 



Monte Carlo [§]| 



b Predicted by conformal invariance |1 



Our results are shown in table 1, where the values of p* and y p for L = 3 and 4 
on the square lattice have been obtained previously []5J. The amplitude of finite-size 
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corrections is much larger than for the corresponding cases of PBC (see e.g. Table I of 
reference IP). However, the finite-lattice sequences are generally well-behaved, allowing 
for smooth extrapolations. Comparing our extrapolated estimates for p c and y p with 
the well-known respective values provides a good overall check of the reliability of our 
procedures. For the square lattice our p c agrees very well with, but is less precise than, 



the best estimate for the percolation threshold p c = 0.592745 ± 0.000002 [[HJ. A similar 
picture holds for the comparison with the exact p c = 1/2 for the triangular lattice, and 
y p = 3/4 (both lattices). 

Turning now to i] s , table 1 provides a direct test of the prediction rj s = 2/3 fi~3H . 
previously confirmed only indirectly via the series result 71 = 2.10 ± 0.02 |l7j which, 
together with the scaling relation 271 = 7 + v(2 — 7] s ) and the exact values 7 = 43/18 
and v — 4/3 gives r] s = 0.64 ± 0.03 . Though the agreement is generally very 

satisfactory, the sequence for the triangular lattice seems to extrapolate towards a region 
slightly above 2/3. 

In order to improve the quality of our estimates |7|, 0, we have also generated 
sequences of finite-size data for r] s by setting p at the best available (or exact) value of 
p c . For comparison with the corresponding data of reference [§] for the bulk exponent 77, 
given by r\ = L/ti^l(p c ) |L2| (where £l is related to the largest eigenvalue of the transfer 
matrix with PBC and, in reference ||, is calculated at the respective p* as in table 1 
above), we have done the same for PBC. The results are displayed in table 2. 

For the square lattice the central estimate of p c from reference |16j has been used. 
Had the respective error bars been taken into account, this would typically give rise to 
uncertainties in the sixth decimal place of rj s or rj. BST extrapolations of the truncated 
values point essentially towards the same limits exhibited in table 2 (though with error 
bars roughly doubled), which shows that the corresponding sequences are rather robust. 

The amplitude of finite-size corrections is much larger for FBC than for PBC, a 
trend already noticed in the discussion of table 1. Comparison with the final estimate 
rj = 0.2088 ± 0.0008 of reference || suggests that, for PBC at least, setting p = p c 
rather than at the approximate p* reduces error bars by one order of magnitude. This 
latter statement assumes that both BST and the extrapolation procedures described in 
reference || are of comparable intrinsic accuracy, which is reasonable for long (~ 10 
elements) sequences such as those encountered here (for shorter sequences, BST would 
be relatively more reliable than other methods |15[)- Turning to FBC, one sees a similar 
improvement in the extrapolated values of r] s relative to those in table 1. Focusing 
on the remarkably smooth sequence for the square lattice, and taking into account the 
effect of the uncertainty in p c on the estimates of i] s as discussed above, we reach the 
final estimate rj s = 0.6664 ± 0.0008. This is entirely in agreement with the prediction 
r] s = 2/3 from conformal invariance |13[], and l| orders of magnitude more accurate 
than previous numerical results [131 |T7|. 
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Table 2. Results for rj s and r\ obtained by setting p — 0.592745 (square) or p = 1/2 
(triangular lattice). Uncertainties in last quoted digits are shown in parentheses. 
Extrapolations obtained by BST algorithm with correction-to-scaling exponent uj in 
ranges shown. Expected values are exact, unless otherwise noted. 

Square Triangular 



L 


Vs 


V 


Vs 


V 


2 


0.4326951882 


0.2163475941 


0.4673843349 


0.2114765825 


3 


0.4848101732 


0.2130595008 


0.5135927770 


0.2111933048 


4 


0.5174925026 


0.2125576128 


0.5423526787 


0.2103549509 


5 


0.5400291833 


0.2114673276 


0.5619831315 


0.2097409643 


6 


0.5565576800 


0.2107370714 


0.5762429521 


0.2093492678 


7 


0.5692194471 


0.2102232886 


0.5870747896 


0.2090956127 


8 


0.5792399531 


0.2098564768 


0.5955841200 


0.2089245599 


9 


0.5873733739 


0.2095868033 


0.6024472809 


0.2088045137 


10 


0.5941102466 


0.2093833099 


0.6081026832 


0.2087173009 


Expected 


2/3 Q 


5/24 


2/3 a 


5/24 


Extrapolated 


0.6664(4) 


0.20835(2) 


0.665(1) 


0.20833(2) 




1.00(5) 


1.90(10) 


1.00(10) 


2.00(5) 



Predicted by conformal invariance [13j 



It is known from the exact solution of the Ising model that finite-size estimates of 
the critical temperature and exponents converge as T C (L) — T c ~ L~ 3 ; y(L) — y ~ L~ 2 
(y = v,r)) for PBC [0, while for FBC the corrections are respectively proportional to 
L~ 2 and L^ 1 |jl9fl . For percolation on strips with PBC, numerical evidence is similarly 
consistent with p c (L) — p c ~ L~ 3 and y p (L) — y p ~ L~ 2 0. In the present case, 
the data of tables 1 and 2 point towards the following scenario: y p {L) — y p ~ L^ 1 , 
T] S (L) —rj s ^ L~ x (FBC); rj(L) — r] ~ L~ 2 (PBC). Though data for the triangular lattice 
indicate p c (L) — p c ~ L~ 2 as expected, for the square lattice one seems to get fits with 
the same quality either for uj = 1 or 2, or just about any value in between. We have 
been unable to sort out this apparently discrepant behaviour. 

We have investigated the possible existence of a higher-order critical point, related 
to a surface-assisted transition. Series work indicates that a special transition should 
not be expected for percolation clusters in two dimensions (though in three-dimensional 
systems it should occur) |T7|, |2T)fl . The work described below is a direct test of such 
results for the two-dimensional case. 

Similarly e.g. to studies of polymer adsorption [lCj, [□]], a distinct occupation 
probability p s is introduced for sites on either strip boundary. Fixed points 
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are obtained by comparing correlation lengths on three strips 0: 

£l(p*,P* s ) = £l-i{p*,P* s ) = £l-2{p*,P* s ) (3) 
L L — l L — 2 ' U 

By analogy with polymer adsorption, if a special transition occurs it must be at some 
p* > p* so that the critical cluster is located predominantly close to the edge. As p* is 
a bulk quantity, one expects it to converge to p c regardless of whether the transition is 
ordinary or special. By scanning the (p,p s ) space we have ascertained that there is only 
one non-trivial solution of Equation (||), which corresponds to an ordinary transition. 
This can be seen from the estimates of critical parameters and respective exponents 
shown in table 3. 

Table 3. Results from two-parameter PR. Uncertainties in last quoted digits are shown 
in parentheses. Extrapolations obtained by BST algorithm with correction-to-scaling 
exponent uj in ranges shown. Expected values are exact, unless otherwise noted. 

(a) Square 



L 


* 

P 


P*s 


Vp 


Vs 


Vs 


5 


0.595339 


0.503680 


0.731343 


- 1.03326 


0.634489 


G 


0.595215 


0.503230 


0.736897 


- 1.04906 


0.635388 


7 


0.594602 


0.500260 


0.740391 


- 1.04301 


0.640568 


8 


0.594148 


0.497414 


0.742624 


- 1.03692 


0.644923 


9 


0.593824 


0.494839 


0.744160 


- 1.03161 


0.648407 


10 


0.593590 


0.492552 


0.745265 


- 1.02719 


0.651168 


Expected 


0.592745(2) a 




3/4 


- 1 


2/3 b 


Extrapolated 


0.5926(1) 


0.460(2) 


0.750(2) 


- 1.001(1) 


0.666(1) 




2.0(4) 


0.85(15) 


2.0(2) 


1.9(1) 


2.00(5) 






(b) Triangular 






L 


P* 


* 

Ps 


Vp 


Vs 


Vs 


5 


0.503487 


0.425365 


0.734458 


- 1.022979 


0.844465 


6 


0.502259 


0.421150 


0.741249 


- 1.027182 


0.802242 


7 


0.501519 


0.417782 


0.744379 


- 1.024663 


0.777288 


8 


0.501064 


0.415120 


0.746081 


- 1.021105 


0.760540 


9 


0.500776 


0.413005 


0.747086 


- 1.017650 


0.748374 


10 


0.500595 


0.411374 


0.747618 


- 1.013932 


0.738967 


Expected 


1/2 




3/4 


- 1 


2/3 b 


Extrapolated 


0.497(2) 


0.402(1) 


0.750(2) 


- 1.004(3) 


0.680(15) 


id 


2.0(4) 


2.0(1) 


2.0(4) 


2.0(4) 


2.0(4) 



Monte Carlo || 



Predicted by conformal invariance |13| 
Once again, the smooth convergence of the sequences of estimates of p*, y p and rj s 
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towards the expected values confirms that our procedures are, on the whole, reliable. 
That the second exponent, y s , is negative ensures that we are dealing with an ordinary 
critical point; our extrapolation is compatible with y s = — 1, a result derived on general 
grounds for the ordinary transition of two-dimensional systems [^TJ . The non-universal 
p* converges to values smaller than the respective p c for each lattice. This resembles the 
ordinary transition for polymers, at which the fugacity for surface contacts is slightly 
de -enhanced pi ITT . 



With the notable exception of the sequence for p* for the square lattice, the leading 
correction-to-scaling exponent seems to be in the vicinity of 2, or even larger, for all 
quantities involved. At present it is not clear whether this feature is fortuitous, or in 
some way related to the structure of the two-parameter PR equations. 

We have shown that the exponent that controls the decay of critical correlations 
along the surface of a semi-infinite percolating plane is rj s = 0.6664 ± 0.0008, consistent 
with the prediction from conformal invariance rj s = 2/3. By setting the site occupation 
probability p at its critical value p c , clean numerical evidence has been provided 
that the finite-size estimates of r] s and of the bulk exponent rj scale respectively as 
f] s {L) — T] s ~ L" 1 (FBC); i](L) — r\ ~ L~ 2 (PBC). It has been shown by numerical 
examination of suitable two-parameter PR recursion relations that no special transition 
occurs in the case; further, at the ordinary critical point the irrelevant exponent is, with 
all probability, y s = — 1 exactly. 

Extensions of the present work to branched polymers (lattice animals) M are 
currently being pursued. Though conformal invariance concepts are not applicable in 
the case P2[ , surface critical indices such as the crossover exponent = y s /y can be 



calculated and compared e.g. to series results [23|, for which error bars are rather large 
at present. 
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